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Abstract 

Photonic events with large missing energy have been observed in e + e _ collisions at centre- 
of-mass energies of 130, 136 and 183 GeV collected in 1997 using the OPAL detector at LEP. 
Results are presented for event topologies with a single photon and missing transverse energy or 
with an acoplanar photon pair. Cross-section measurements are performed within the kinematic 
acceptance of each selection. These results are compared with the expectations from the Standard 
Model process e + e _ — * vv + photon(s). No evidence is observed for new physics contributions to 
these final states. Using the data at ^/s = 183 GeV, upper limits on cr(e + e~ — > XY) • BR(X — » Y7) 
and cr(e+e- -> XX) ■ BR 2 (X -> Y7) are derived for the case of stable and invisible Y. These limits 
apply to single and pair production of excited neutrinos (X = v* , Y = u), to neutralino production 
(X = X2, Y = Xi) and to supersymmetric models in which X = Xi and Y = G is a light gravitino. 
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1 Introduction 



We describe measurements and searches using a data sample of photonic events with large missing 
energy collected in 1997 with the OPAL detector at LEP. The events result from e + e~ collisions 
at centre-of-mass energies of 130.0, 136.0 and 182.7 GeV with integrated luminosities of 2.35, 3.37 
and 54.5 pb -1 , respectively. The present paper complements our recent publication of results from 
earlier data samples Q consisting of a total of 25 pb -1 at centre-of-mass energies of 130, 136, 161 
and 172 GeV. Results on photonic events without missing energy at ^/s = 183 GeV are presented in 
a separate paper |2| . Recent measurements of photonic event production have also been made by the 
other LEP collaborations at centre-of-mass energies above the W pair threshold Q, including new 
results at yfs = 183 GeV [||. Related searches in pp collisions have been reported in H. The interest 
in the 1997 data is twofold. The main motivation is that the large data set at a higher centre-of-mass 
energy (183 GeV) gives discovery potential in a new kinematic regime. Additionally, the lower energy 
data sets allow the experiment at 130 and 136 GeV to be repeated leading to a check of our previous 
results which indicated an excess of events with acoplanar photons 

The single-photon and acoplanar-photons search topologies presented here are designed to select 
events with one or more photons and significant missing transverse energy, indicating the presence 
of at least one neutrino- like invisible particle which interacts only weakly with matter. The event 
selections for these search topologies are similar to those used in our recent publication Q. The main 
changes are improvements in the rejection of backgrounds from cosmic ray interactions which allow for 
an increased kinematic acceptance of the single-photon topology. The single-photon search topology 
is sensitive to neutral events in which there are one or two photons and missing energy which, within 
the Standard Model, are expected from the e + e~ — > 1/1/7(7) process^. The allowance for two observed 
photons retains acceptance for doubly radiative neutrino pair production. The acoplanar-photons 
search topology is designed to select neutral events with two or more photons and significant missing 
transverse energy which, within the Standard Model, are expected from the e + e~ — ► 1/1/77(7) process. 
The selection is designed to retain acceptance for events with three photons, if the system formed by 
the three photons shows evidence for significant missing transverse energy. 

These photonic final-state topologies are sensitive to new physics of the type e + e~ — * XY and 
e + e~ — ► XX where X is neutral and decays radiatively (X — > Y7) and Y is stable and only weakly 
interacting. For the general case of massive X and Y this includes conventional supersymmetric 
processes || (X = X2, Y = x?)- In this context it has been emphasized that the radiative branching 
ratio of the %\ m ay be large. These topologies also have particularly good sensitivity for the special 
case of My ~ 0. This applies both to the production of excited neutrinos (X = u*,Y = v) and to 
supersymmetric models in which the lightest supersymmetric particle (LSP) is a light gravitincQ and 
Xi is the next-to-lightest supersymmetric particle (NLSP) which decays to a gravitino and a photon 
(X = xJ,Y = G). For supersymmetric models with a light gravitino, the photonic branching ratio of 
the xi is naturally large. Such a signature has been discussed in || and more recently in ||-12| in the 



context of both no-scale supergravity models and gauge-mediated supersymmetry breaking models. 
Other types of new physics to which these search topologies are sensitive include the production of 
an invisible particle in association with a photon and the production of invisible particles tagged 
by initial-state radiation. An example of this is e + e~ — > GG7 |13}| . The acoplanar-photons search 
topology also has sensitivity to the production of two particles, one invisible, or with an invisible decay 
mode, and the other decaying into two photons. Such events might arise from the production of a 
Higgs-like scalar particle, S° : e + e~ — ► Z°S°, followed by S ^ 77, Z° — > vu. OPAL results from a 
search for this process at = 183 GeV, including the hadronic and leptonic Z° decays, are reported 



1 The photon in parentheses denotes that the presence of this photon is allowed but not required. 
2 The mass scale is typically O(keV). 
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in a separate paper |fl4j1 . The results from our previous searches |IJ for e + e _ — > XY and e + e _ — > XX 
with X —* Y7 have been used to set model-dependent limits on excited neutrinos and neutralinos [jig ]. 
The new results reported here can be used in the same manner. 

This paper will first briefly describe the detector and the Monte Carlo samples used. The event 
selection for each search topology will then be described, followed by cross-section measurements for 
e + e~ — > 1^7(7) and e + e~ — > ^1777(7) and comparisons with Standard Model expectations. Implica- 
tions of these results on the possibility of new physics processes will be discussed. 



2 Detector and Monte Carlo Samples 



The OPAL detector, which is described in detail in [16], contains a pressurized central tracking system 
operating inside a solenoid with a magnetic field of 0.435 T. The region outside the solenoid (barrel) and 
the pressure bell (endcap) is instrumented with scintillation counters, presamplers and the lead-glass 
electromagnetic calorimeter (ECAL). The magnet return yoke is instrumented for hadron calorimetry 
and is surrounded by external muon chambers. Electromagnetic calorimeters close to the beam axis 
measure luminosity and complete the acceptance. 

The measurements presented here are mainly based on the observation of clusters of energy de- 
posited in the lead-glass electromagnetic calorimeter. This consists of an array of 9,440 lead-glass 
blocks in the barrel (| cos#| < 0.82) with a quasi-pointing geometry and two dome-shaped endcap ar- 
rays, each of 1,132 lead-glass blocks, covering the polar angleQ range (0.81 < | cos#| < 0.984). Energies 
measured in the ECAL will normally refer to those obtained after corrections to account for losses in 
upstream material. Energies without these corrections are called deposited energies. In the regions 
(0.72 < I cos# I < 0.82) and (| cos 6 \ > 0.945), the energy resolution of the ECAL is degraded relative 
to the nominal resolution. This degradation is largely due to increased amounts of material in front 
of the ECAL. In some cases (where stated) these regions have been excluded from the analysis. Fully 
hermetic electromagnetic calorimeter coverage is achieved beyond the end of the ECAL down to small 
polar angles with the use of the the gamma-catcher calorimeter, the forward calorimeter (FD) and 
the silicon-tungsten calorimeter (SW). These detectors cover the angular regions with respect to the 
beam of 140-205 mrad, 40-145 mrad and 24-59 mrad, respectively. However, a small region centred on 
a polar angle of 30 mrad lacks useful calorimetric coverage due to the installation, in 1996, of a thick 
tungsten shield designed to protect the tracking chambers from synchrotron radiation background. 
The effective limit of electromagnetic hermeticity is therefore around 33 mrad. 

Scintillators in the barrel and endcap regions are used to reject backgrounds from cosmic ray 
interactions and to provide time measurements for the large fraction (~ 80%) of photons which convert 
in the material in front of the ECAL. The barrel time-of-flight (TOF) scintillator bars are located 
outside the solenoid in front of the barrel ECAL and match its geometrical acceptance (| cos 9\ < 0.82). 
Arrays of thin scintillating tiles with embedded wavelength-shifting-fibre readout [[l?]] were installed 
prior to the 1997 run. The new tile endcap (TE) scintillator arrays are located at 0.81 < | cos 6 \ < 0.955 
behind the pressure bell and in front of the endcap ECAL. 

Additional scintillating tile arrays, referred to as the MIP-PLUG, were installed at polar angles 
between 40 and 200 mrad and consist of two pairs of scintillating tile layers designed for detection of 
muons. The outer pair, covering the angular range of 125 to 200 mrad is used in this paper to provide 



3 In the OPAL coordinate system, is the polar angle defined with respect to the electron beam direction and 
the azimuthal angle. 
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redundancy in the rejection of events with energetic electromagnetic showers in the gamma-catcher 
region. 



The tracking system, consisting of a silicon microvertex detector (SI), a vertex drift chamber (CV) 
and a large volume jet drift chamber (CJ), is used to reject events with prompt charged particles. 
The silicon microvertex detector consists of two concentric cylindrical layers of silicon microstrip 
arrays, each layer providing both an azimuthal and longitudinal (along the beam direction) coordinate 
measurement. The two-layer acceptance covers | cos#| < 0.90. 

Beam related backgrounds and backgrounds arising from cosmic ray interactions are rejected using 
the scintillator timing measurements and information from the electromagnetic calorimeter shower 
shape, the hadron calorimeter and the muon detectors. The integrated luminosities of the data samples 
are determined to better than 1% from small-angle Bhabha scattering events in the SW calorimeter. 



Triggers [18| based on electromagnetic energy deposits in either the barrel or endcap electromagnetic 
calorimeters, and also on a coincidence of energy in the barrel electromagnetic calorimeter and a hit in 
the TOF system, lead to full trigger efficiency for photonic events passing the event selection criteria 
described in the following section. 

For the expected Standard Model signal process, e+e" — * vV + photon(s), the Monte Carlo 
generator KORALZ pll was used. For other expected Standard Model processes, a number of different 



generators were used: RADCOR |§ for e+e" -> 77(7); BHWIDE |2| and TEEGG [|| for e+e 



e + e"(7); grc4f pi for e+e" -> £+£"^1/(7); and KORALZ for e+e" -> //+//" (7) and e+e" -> t + t~(j). 



The expected contributions from each of these Standard Model processes were evaluated using a total 
equivalent integrated luminosity at least ten times larger than the integrated luminosity of the data 
sample. 

To simulate possible new physics processes of the type e+e" — ► XY and e+e" — > XX where X decays 
to Y7 and Y escapes detection, a modified version of the SUSYGEN [^J] Monte Carlo generator was 
used to produce neutralino pair events of the type e+e" — > X2X1 an d e+e" — > X2X2' x\ X?7> with 
isotropic angular distributions for the production and decay of x\ an d including initial-state radiation. 
Monte Carlo events were generated at 40 (for XY production) and 42 (for XX production) points in 
the kinematically accessible region of the (Mx, My) plane. All the Monte Carlo samples described 
above were processed through the OPAL detector simulation p5| . 



3 Photonic Event Selection 



This section describes the criteria for selecting single-photon and acoplanar-photons events. The 
kinematic acceptance of each selection is defined in terms of the photon energy, £L, and the photon 
polar angle, 8. In addition, the scaled energy, x 7 , is defined as Ey/Eb eam , and the scaled transverse 
energy, xt, as x 7 sin#. 



Single-Photon - One or two photons accompanied by invisible particle(s): 

• At least one photon with x T > 0.05 and with 15° < 9 < 165° (| cos 6* | < 0.966). 
Acoplanar-Photons - Two or more photons accompanied by invisible particle(s): 

• At least two photons, each with x 7 > 0.05 and 15° < < 165°, or one photon with JSL > 
1.75 GeV and I cos 6 \ < 0.8 and a second photon with £ 7 > 1.75 GeV and 15° < 9 < 165°. 
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• The scaled transverse momentum of the two-photon system consisting of the two highest 
energy photons, p~f ( , must satisfy pJV-Ebeam > 0.05. 

In each of the two cases, it is desirable to retain acceptance for events with additional photons, 
if the resulting photonic system is still consistent with the presence of significant missing energy. 
This reduces the sensitivity of each measurement to the modelling of higher-order contributions. 
Consequently, a large fraction of the kinematic acceptance of the acoplanar-photons selection is also 
contained in the kinematic acceptance of the single-photon selection. 

3.1 Single-Photon Event Selection 

The single-photon selection criteria are similar to the previous OPAL analysis of photonic events with 
missing energy [jlj but have increased acceptance for lower energy photons due to improved rejection 
and control of cosmic ray and beam related backgrounds. The modifications also improve the efficiency 
for events with two detected photons. 

• Kinematic acceptance. Events must contain a primary electromagnetic cluster (that with 
the highest deposited energy in the barrel or endcap calorimeters) in the region 15° < 9 < 165° 
(| cos 9 | < 0.966) with xt > 0.05. Events are considered to have more than one photon if 
additional electromagnetic clusters are found in the barrel or endcap calorimeter (|cos0| < 
0.984) having deposited energy exceeding 300 MeV. 

• Cluster quality. The primary electromagnetic cluster, combined with any clusters contiguous 
with it, must be consistent with the cluster size and energy sharing of blocks for a photon 
coming from near the interaction point. The cluster size varies in both azimuthal and polar angle 
extent as a function of | cos 9 \ . The cluster extent cuts are parametrized in | cos 9 | accordingly. 
Events are rejected if the cluster energy exceeds the beam energy by more than three standard 
deviations. 

• Forward energy vetoes. Events are rejected if the energy sum of gamma-catcher clusters in 
either end is greater than 5 GeV. Events are also rejected if either layer of the outer MIP-PLUG 
shows evidence of an energetic shower (pulse-height exceeding about ten minimum ionizing 
particle equivalents). An event can also be vetoed based on the transverse momentum sum 
of clusters measured in the forward calorimeters FD and SW. Events are rejected that have a 
transverse momentum sum exceeding 1 GeV and where the azimuthal angle of the transverse 
momentum sum is within 60° in azimuth of the direction opposite the measured momentum 
of the photonic system. A final complementary veto rejects events that have an energy sum 
exceeding 5 GeV where the sum is over all clusters in the FD, SW and forward part of the 
endcap ECAL (| cos 9 \ > 0.966) which are within 60° in azimuth of the direction opposite the 
measured momentum of the photonic system. The directional nature of these last two vetoes 
removes events with forward going high-energy particles that can account for some or all of the 
missing transverse momentum. However, it minimizes losses from random noise or accidental 
energy deposits in the forward detectors. It also allows for the presence of initial-state radiated 
photons in the forward region which are not back-to-back with the photon(s). 

• Muon veto. Events are rejected if there are any muon track segments reconstructed in the 
muon chambers or in the hadron calorimeters. Events are also rejected if there is significant 
activity in the outer part of the barrel hadron calorimeter. The muon veto is used primarily to 
remove cosmic ray background. 
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• Selective multi-photon veto. This veto addresses backgrounds, principally from e + e~ — ► 
77(7)5 whilst retaining acceptance for events with two photons and missing energy. Events with 
a second photon are rejected if any of the following criteria are satisfied: 

— The total energy of the two clusters exceeds 0.9-^5. 

— The acoplanarity anglef] of the two clusters is less than 2.5°. 

— The missing momentum vector calculated from the two clusters satisfies | cos# m i ss | > 0.9. 

— A third electromagnetic cluster is detected with deposited energy exceeding 300 MeV. 

— pV/E heaia < 0.05. 

— For events with at least one of the two clusters in the region | cos 9 \ > 0.95, the variable &t 
is less than 0.1, where bx = (sin#i + sin #2)! cos [(^i — 4>2)/'A |- This amounts to a stronger 
acoplanarity cut for events with at least one forward photon. 

— Events with an energy sum greater than 1 GeV in either end of FD or SW are rejected 
if the two-photon plus forward photon system is planar, namely if the sum of the three 
opening angles exceeds 350°. 

Events are required to either contain no reconstructed charged tracks or contain a photon candidate 
consistent with a photon conversion observed within the central tracking volume. These events are 
referred to as non-conversion and conversion candidates, respectively. These two classes are mutually 
exclusive. We will now describe in turn the additional criteria used for each candidate class. 

Non-conversion candidates must satisfy the following additional criteria : 

• Charged track veto. It is required that there are no reconstructed tracks with 10 or more hits 
in CJ. 

• Timing requirements. The criteria depend on the polar angle of the primary cluster. For the 
angular region | cos# | < 0.72, either an in-time associated TOF hit or the absence of the special 
background vetoes (described below) is required. The second requirement is made in order to 
retain acceptance for photons that do not convert before reaching the TOF counters. For the 
range 0.72 < | cos9 | < 0.82, where the muon coverage is not complete and the material in front 
of the ECAL leads to a high probability of a photon being detected in the TOF, a good in-time 
associated TOF hit is required. For | cos6> | > 0.82, a good in-time associated TE hit is required 
as well as the absence of the first three of the special background vetoes described below. A 
cluster with an associated TOF (TE) hit is considered to be in-time if the measured arrival time 
of the photon at the TOF (TE) is within 5 (30) ns of the expected time for a photon originating 
from the interaction point. These definitions of in-time hits also apply to the single-photon 
conversion selection as well as the acoplanar-photons selection. Events with a photon candidate 
having an out-of-time associated TOF or TE hit are rejected as cosmic rays. 

• Special background vetoes. Three special background vetoes are used for candidates that are 
in the endcap region or that have no TOF timing information. A fourth special background veto 
is used only for photon candidates with | cos 9 \ < 0.72 that have no TOF timing information. The 
first veto rejects events in which any of the three muon triggers (barrel and two endcaps) were 
present. This veto rejects cosmic ray background. The second looks for a series of electromagnetic 
or hadronic calorimeter clusters consistent with the same radial and azimuthal position as the 
primary cluster, but at different positions along the beam direction. This veto rejects beam halo 

4 Defined as 180° minus the opening angle in the transverse plane. 
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backgrounds. The third looks for a series of hits in the outer layers of the hadron calorimeter. 
This veto rejects both cosmic rays and beam related backgrounds. The fourth veto is based 
on the shape of the cluster in the barrel ECAL. The observed energy deposition in each lead- 
glass block of the cluster is fitted to the expected shower profile for a photon coming from 
the interaction point. One then calculates the variable S (in radiation lengths), defined as the 
difference between the measured and expected values of the energy-weighted lateral distance 
from the fitted shower centroid. Events with S exceeding 0.2 are rejected. This veto rejects 
cosmic rays and beam related backgrounds, both of which tend to have shower shapes wider 
than those of photons originating from the interaction point. 

Conversion candidates must satisfy the following additional criteria : 

• Photon conversion consistency. There must be at least one reconstructed charged track 
in the central tracking chambers. The charged track with the most hits must be associated in 
space to at least one of the two most energetic photon candidates. In particular, the measured 
polar angle of the track should be consistent with the polar angle of the ECAL cluster to within 
100 mrad, and the azimuth at the point of closest approach of the track to the interaction point 
should differ by less than 100 mrad from the measured azimuthal angle of the ECAL cluster. 
The azimuthal matching criterion is relaxed to 500 mrad for clusters with |cos#| > 0.90 since 
conversions at forward angles often lead to large showers and difficulties in resolving the jet 
chamber left-right ambiguity. In addition, for such forward photon conversion candidates, which 
mostly convert in the CV endplate, it was further required that there be at least 12 out of the 
first 16 wires hit in CJ. 

• Prompt charged track veto. Events are rejected as being consistent with containing a prompt 
charged track if at least one photon candidate has azimuthally associated hits in the innermost 
tracking detector (SI for | cos# | < 0.9 and CV otherwise). 

• Two or more track veto. Events with conversion candidates are rejected if they have at least 
two tracks, reconstructed from axial-wire hits in CV, with an opening angle in the transverse 
plane exceeding 45°. This criterion is used principally for redundancy in the rejection of Bhabha 
scattering events. 

• Identified cosmic ray veto. Events with at least two electromagnetic clusters in the barrel 
region each with associated TOF hits are rejected as identified cosmic rays if the time difference 
between the upper and lower TOF hits is consistent with a downward-going cosmic ray. 

• Timing requirements. The photon(s) associated with the track is required to have an in-time 
TOF or TE hit depending on whether the polar angle of the photon matches the TOF or TE 
geometrical acceptance. 

Distributions of some quantities used in the event selection are illustrated in Figures and |. 
For photons with | cos 8 \ < 0.82, Figure [l|a shows the difference between the measured TOF timing 
and that expected for a photon originating from the interaction point for events passing all selection 
criteria or failing only the TOF timing requirement. The eleven events outside the accepted region of 
±5 ns are rejected as cosmic rays. The expected region for good events is shown in greater detail in 
Figure [l|b. Figure [l|c shows the corresponding plot for photons measured in TE with | cos 6 \ > 0.82. 
The three events outside the accepted region of ±30 ns are rejected as cosmic rays. Figure |2|a shows 
the distribution of the cluster shape variable S for events passing all selection criteria or failing only the 
cluster shape cut of the special background vetoes. Rejecting events with S greater than 0.2 preserves 
the e + e~ — > 1/177(7) Monte Carlo events, while in the data it removes very large clusters which have 
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been verified as being due to cosmic rays and beam related backgrounds. Figure [|b shows the effect 
of one of the directional forward veto cuts against its principal intended background, e + e~ — > e + e~j. 
This cut is designed to remove events in which a highly energetic particle travelling close to the beam 
direction balances most or all of the transverse momentum of the observed photon(s). The events 
plotted are those passing all cuts or failing only this cut. All the remaining e + e _ — ► e + e~(7) Monte 
Carlo background events fall into this cut region while the signal from the e + e~ — > 1^7(7) Monte 
Carlo has only a very small number there. The data distribution for this plot was checked and found 
to agree well with the Monte Carlo prediction. 



3.2 Acoplanar-Photons Event Selection 

The acoplanar-photons selection has two overlapping regions of kinematic acceptance in order to retain 
both sensitivity to low-energy photons and acceptance at large | cos 9 \ . These selections are based on 
analyses previously published by OPAL using data collected at centre-of-mass energies of 130-172 
GeV Q. Except where specified, cuts on photon candidates apply to the two highest energy photon 
candidates found within the kinematic acceptance. The event selection criteria are described below: 



• Kinematic acceptance. Events are accepted as candidates if there are at least two electro- 
magnetic clusters with scaled energy, x 7 , exceeding 0.05 in the polar angle region 15° < 9 < 165° 
(| cos 9 I < 0.966). In order to retain sensitivity to physics processes producing low-energy pho- 
tons, the minimum energy requirement is relaxed to 1.5 GeV deposited energy (corresponding 
to a photon energy of about 1.75 GeV [26]) for events with a photon candidate in the polar 



angle region | cos 9 \ < 0.8 provided this photon is associated to an in-time TOF hit as outlined 
below in the description of the timing requirements. These two selections are referred to below 
as the "high-energy" and "low-energy" selections, respectively. Background vetoes are applied 
differently for the two parts of the selection, as described below. The system consisting of the 
two highest energy photons must satisfy p^ 7 /-Ebcam > 0.05. 

• Photon conversion consistency requirements or charged track veto. For the high- 
energy selection, events having tracking information consistent with the presence of at least 
one charged particle originating from the interaction point are rejected. The rejection criteria 
are designed to retain acceptance for events in which one or both of the photons convert. Hit 
information from each of CJ, CV, and SI (for | cos 9 \ < 0.9) are used to form independent 
estimators for the existence of charged particle activity. Events are rejected on the basis of 
azimuthal association of charged particle activity with the photon candidate clusters. To reject 
e + e~ — ► £ + ^~77, an additional veto requires that there be no reconstructed charged track with 
transverse momentum exceeding 1 GeV, with associated hits in CV, and separated from each of 
the photon candidates by more than 15°. 

The low-energy part of the selection does not allow photon conversions in the tracking chambers. 
It requires that there be no reconstructed charged track in the event with 20 or more hits in CJ. 

• Cluster quality. Photon candidates within the polar angle region | cos 9 | < 0.75 are required 
to have an angular cluster extent less than 250 mrad in both 9 and <fi. Additionally, to reduce 
background from cosmic rays which graze the electromagnetic calorimeter producing extended 
energy deposits that can be split by the clustering algorithm, photon candidates are required 
to be separated by at least 2.5° in azimuth. Events are rejected if a photon candidate cluster 
energy exceeds the beam energy by more than three standard deviations. 

• Forward energy vetoes. The forward vetoes described for the single-photon selection are 
applied with the same thresholds. 
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• Muon veto. To suppress backgrounds arising from cosmic ray muon interactions or beam halo 
muons which can deposit significant energy in the calorimeter, the events must pass the muon 
veto described for the single-photon selection. Additionally, the first three individual vetoes of 
the special background vetoes described for the single-photon selection are applied to events in 
which no TOF information is present. 

• Timing requirements. For the low-energy part of the selection, in order to ensure that the 
trigger is fully efficient for low-energy photons, we require that there be a photon in the barrel 
region with an associated in-time TOF hit. For the high-energy selection the event must have an 
associated in-time TOF or an associated in-time TE hit for at least one of the photon candidates. 
Events with a photon having an associated out-of-time TOF hit are rejected. Events in which 
one photon candidate has an associated out-of-time TE hit are retained provided they pass the 
other timing requirements. Finally, if there is a charged track associated with a cluster within the 
polar angle region | cos 9 \ < 0.82, the requirement of an associated in-time TOF hit is applied. 

• Selective multi-photon veto. As for the single-photon selection, this veto is designed to 
reject backgrounds primarily from e + e~ — ► 77(7) whilst retaining acceptance for events with 
two or more photons and missing energy. Events are rejected if any of the following criteria are 
satisfied: 

— The total visible ECAL energy of the event exceeds 0.95-^/5. 

— The acoplanarity angle of the two highest energy clusters is less than 2.5°. 

— The missing momentum vector calculated from the two highest energy photon candidates 
satisfies I cos m i ss I > 0.95. 

— Events having three or more photon candidates (with deposited energy greater than 300 MeV) 
are rejected unless the system formed by the three highest energy photons is significantly 
aplanar (sum of the three opening angles < 350°) and the transverse momentum of the 
three-photon system exceeds O.lE'beam- F° r events with an energy sum greater than 1 GeV 
in either end of FD or SW, the aplanarity cut is applied, using the forward detector as the 
third photon candidate. 



The acoplanar-photons selection described above has a lower energy threshold for the most ener- 
getic photon than the single-photon selection. However the single-photon selection has more accep- 
tance for events without time-of-flight information for the photons. In order to obtain the best overall 
acceptance for acoplanar-photons, we have added to the above described acoplanar-photons selection 
that part of the single-photon selection which is within the kinematic acceptance of the acoplanar- 
photons selection. This addition results in a relative increase in efficiency of 9% for Standard Model 
e + e~ — > 1/1/77(7) events. 



4 Results 



The results of the single-photon and acoplanar-photons selections are given below in sections iA 
and 4.2. The measured cross-sections for each search topology are given and compared with Standard 
Model expectations. As no evidence for new physics processes is seen, the results are presented in 
terms of upper limits on a(e + e~ XY) • BR(X -> Y7) and <7(e+e _ XX) • BR 2 (X -> Y7). This is 
done both for the general case of massive X and Y, applicable to conventional super symmetric models 
in which X = X2 an d Y = %?, and also separately for the special case of My ~ 0, which applies both 
to single and pair production of neutralinos in supersymmetric models in which the LSP is a light 
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gravitino and to single and pair production of excited neutrinos. All efficiencies are evaluated with 
the decay length of X set to zero. 

For the purposes of new physics searches, only the y/s = 183 GeV data are considered; the 
i/s = 130 and 136 GeV data do not open any new kinematic regions, nor is the integrated luminosity 
at these energies sufficient to significantly improve the potential for discovery. For both the XX and 
XY searches, Monte Carlo samples were generated for a variety of mass points in the kinematically 
accessible region of the (Mx, My) plane. To set limits for arbitrary Mx and My, the efficiency over the 
entire (Mx,My) plane is parametrized using the efficiencies calculated at the generated mass points. 
In the single-photon search topology, the regions with Mx + My < M% are kinematically accessible 
at y/s ~ Mz, and strong limits have already been reported pTfl . In the acoplanar-photons search 
topology, limits have been reported for masses Mx < Mz/2 Pq| . In these low mass regions, possible 
radiative return to the Z° followed by Z° — » XY or XX would yield very different event kinematics than 
those produced by the signal Monte Carlo generator. For these reasons, the search for XY production 
is restricted to the region with Mx + My > Mz, and the search for XX production is restricted to 
Mx values larger than about Mz/2. 

4.1 Single-Photon 

After applying the selection criteria of the single-photon selection to the yfs = 130, 136 and 183 GeV 
data samples, a total of 21, 39 and 191 events are selected. The expected contributions from cosmic 
ray and beam related backgrounds are 0.02, 0.02 and 0.4 events, respectively. These backgrounds have 
been estimated from events having out-of-time TOF or TE information but passing all other selection 
criteria and from events selected with looser criteria that have been visually scanned. Of the expected 
physics backgrounds from plausible sources, only e + e~ — > 77(7), e + e~ — > £ + l~vD(^), e + e~ — > [i + /i~7 
and e + e~ — > r + r~7 have non-negligible contributions. Coming primarily from the e + e~ — ► ~i~ 
and fj, + n~'y final states, the total physics backgrounds contribute 0.07, 0.09 and 0.4 events, respectively, 
to the 130, 136 and 183 GeV samples. The background contributions are summarized in Table [I]. For 
each of the three centre-of-mass energies, Table § shows the number of events observed, the number 
of events expected from the Standard Model process e + e~ — > 1/1/7(7) evaluated using the KORALZ 
generator and the sum of background events expected from other Standard Model physics processes 
with those from cosmic ray and beam related processes. The numbers of events observed agree with 
the numbers expected from e + e~ — > 1/7/7(7) plus the background. The estimated efficiencies for 
selecting e + e~ — > 1/1/7(7) events within the kinematic acceptance of the single-photon selection are 
also given in Table |2|, as are the corresponding measured e + e~ — > 1/1/7(7) cross-sections within this 
kinematic acceptance, corrected for detector and selection efficiencies, and subtracting the estimated 
background. For both the single-photon and acoplanar-photons selections, efficiency losses due to 
vetoes on random detector occupancy range from about (2-4)% at the different centre-of-mass energies. 
Quoted efficiencies include these losses. 

The total systematic error on the cross-section measurement is estimated to be 3.5%. The con- 
tributing uncertainties are from the integrated luminosity (0.5%), effects due to uncertainties on the 
electromagnetic calorimeter energy scale and resolution (0.7%) and the detector occupancy estimate 
(1%). In addition, the overall selection efficiency uncertainty (2%) is caused mainly by uncertainties in 
the simulation of the detector material and consequent photon conversion probabilities. An additional 
error of 2.5% is assigned based on the comparison of the estimated efficiency using two different event 
generators [fl]]. The cross-section as a function of centre-of-mass energy is plotted in Figure ^. In 
this plot the measured cross-sections for yfs = 130 and 136 GeV are 11.1 ± 1.7 and 15.9 ± 1.9 pb, 
respectively. These are weighted averages of the 1995 data results with the results from this analysis. 
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The results for y/s = 161 GeV and \fs = 172 GeV are also plotted. The cross-section results from 
our earlier publication |l|] have been corrected for the slightly different kinematic acceptance used in 
the analysis of those data. The corrected cross-sections, obtained for the 1995 and 1996 data at y/s = 
130, 136, 161 and 172 GeV are 10.6 ± 2.4, 17.3 ± 3.0, 5.6 ± 0.8 and 5.8 ± 0.8 pb, respectively. For 
1997 data at \fs = 130, 136 and 183 GeV the measured cross-sections are 11.6 db 2.5, 14.9 db 2.4 and 
4.71 d= 0.34 pb, respectively. The curve shows the predicted cross-section from the KORALZ event 
generator for the Standard Model process e + e~ — * 1/1/7(7). The data are in reasonable agreement with 
the prediction. 

In Figure ||a, the scaled energy of the most energetic photon is plotted against the cosine of its 
polar angle for events in the ^/s=183 GeV sample. The data are distributed as expected from the 
e + e~ — > 1/1/7(7) Monte Carlo. Similar agreement is found for the 130 and 136 GeV data. In Figure ||b 
the polar angle distribution for the ^/s=183 GeV sample is shown and agrees with the e + e~ — * 1/1/7(7) 
Monte Carlo expectation. If one calculates the recoil mass M reco n, defined as the mass recoiling against 
the photon (or against the two-photon system), one expects a peak in the M recoi i distribution at Mz, 
due to a large contribution from the decay Z° — ► uu. One clearly sees this feature in the data as shown 
in Figure There is good agreement between data and Monte Carlo in this distribution for each of 
the three centre-of-mass energies. 

The single-photon selection was designed to allow for the presence of a second photon in order to 
accept events from the e + e~ — ► z/I/77 process. In the \/s = 130, 136 and 183 GeV data sets, 2, 4 and 
12 observed events are considered to be two-photon events (i.e. have a second photon with deposited 
energy exceeding 300 MeV in the ECAL). This is consistent with the expectations of 1.4, 1.7 and 11.3 
events, respectively, from the KORALZ Monte Carlo. 

4.1.1 Search for e+e - — > XY, X — > Y7 ; General case: My > 

The single-photon selection described in Section ^ is designed to maximize acceptance for Standard 
Model e + e~ — > 1/1/7(7) events. However, when searching for signatures of the process e + e~ — ► XY, X — ► 
Y7, it is possible to implement further cuts to reduce the contribution from e + e~ — > 1/1/7(7), which is 
now considered as background. Depending on the values of Mx and My, various combinations of the 
following cuts are applied to events in the single-photon sample: 

• Kinematic consistency: The energy of the most energetic photon is required to lie within 
the range kinematically accessible to a photon from the process e + e~ — * XY, X — ► Y7, after 
accounting for energy resolution effects. 

• Degraded resolution: Events with M reco ji <75 GeV are rejected if the most energetic photon 
lies in the region 0.72 < | cos 6* | < 0.82 or |cos#| > 0.945. Energy resolution in these angular 
regions is known to be degraded. 

• Low recoil mass: Require M reco ii <75 GeV. 

• Z° radiative return: The event is rejected if 75 GeV < M recoi i < 105 GeV. 

The kinematic consistency cut and the degraded resolution cut are applied to all (Mx, My) values. 
The low recoil mass cut is applied if My < 0.28Mx — 18 (GeV), in which case a significant portion of 
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the expected photon energy distribution is consistent with low recoil masses^]. The radiative return 
cut is applied if the low recoil mass cut is not applied, and if Mx and My are such that the difference 
between the maximum and minimum kinematically allowed photon energies (before energy resolution 
effects are considered) is greater than 0.3i?beam- 

The selection efficiencies including the above cuts as a function of (Mx, My) are given in Figure ||. 
Also illustrated is the region in which the radiative return cut is applied. In the region where Mx — My 
is small, photon energies are correspondingly small and efficiencies become low. Since uncertainties 
due to energy scale and resolution effects become significant in this region, we do not consider values 
of Mx and My that lead to efficiencies of less than 40% in the absence of the recoil mass or radiative 
return cuts. 

The number of selected events in the data consistent with each (Mx, My) value is shown in Figure |7| 
and can be compared with the number expected from Standard Model e + e~ — > 1/1/7(7) events as shown 
in Figure ||. In general there is good agreement, and we proceed to set upper limits at 95% confidence 
level (CL) on the cross-section times branching ratio, <r(e + e - — > XY) • BR(X — > Y7), which are shown 
in Figure |9[ The upper limits are calculated taking into account the expected number of Standard 
Model e + e~ — > 1/1/7(7) background events estimated from KORALZ using the method described 
in p9|| . Background from sources other than e + e~ — ► 1/1/7(7), including the estimated cosmic ray 
and beam related background, is intentionally not taken into account in the limit calculations. The 
resulting upper limits range from 0.075 pb to 0.80 pb. 

The systematic error on the efficiency for selecting events from potential new physics sources is due 
to the effects already discussed in section 4.1, as well as the uncertainty on the efficiency parametriza- 
tion across the (Mx,My) plane. The parametrization was compared to efficiencies obtained from 
the fully simulated Monte Carlo samples at 40 selected (Mx,My) points, and a resulting systematic 
error of 1% (absolute) was assigned across the plane. The total relative systematic error varies from 
2 to 6% depending on Mx,My; its effect on the upper limits is small, and is calculated according 



to [30]. Uncertainties in the 1/1/7(7) background estimate are also taken into account. Unlike the 
systematic error on the efficiency, a relatively small uncertainty in the background estimate can have 
a significant impact on the resulting 95% CL limit, especially when the number of expected events is 
large. To account for this uncertainty, a convolution is performed within the upper limit calculations. 
Contributing sources to the background uncertainty are: the factors considered in the cross-section 
measurement (3.5%), the estimated theoretical uncertainty in the 1/1/7(7) cross-section (2%) based in 
part on a comparison between the KORALZ and NUNUGPV98 |Il]] event generators, and a contri- 
bution due to uncertainties in the energy scale. This last contribution is dependent on the values of 
Mx and My; it is calculated separately at each (Mx, My) point, and ranges from negligible to 5%. 



4.1.2 Search for e+e - — > XY, X — > Y7 ; Special case: My « 

The case My ~ is applicable to excited neutrino models and to some supersymmetric models 
mentioned earlier. The results presented above include this case and no separate analysis is performed, 
but the results are highlighted here. As described earlier, the 75 GeV recoil mass cut is applied 
for all Mx hypotheses, so that the expected number of events is small. For example, in the range 
91 < Mx < 170 GeV, the expected contribution from 1/1/7(7) is 0.98 ± 0.12 events and there is one 
event observed. Although the numbers of expected and observed events are constant in this range, the 

5 This region was chosen so as to optimize the expected sensitivity. The optimization condition chosen was that the 
expected upper limit on a ■ BR for new physics contributions be minimized, where the expected upper limit is defined as 
the average limit one would expect to set in the absence of new physics contributions. This definition has the advantage 
that it does not require one to specify the cross-section of possible new physics. 



13 



efficiency of the recoil mass cut increases with increasing values of Mx, leading to decreasing values 
for the resulting upper limits on cr(e + e~ — * XY) • BR(X — ► Y7). For Mx > 170 GeV, the kinematic 
consistency requirements become more restrictive than the recoil mass cut. There are no longer any 
events kinematically consistent with hypotheses of Mx > 173 GeV; the background expectation in 
this region varies from 0.60 ±0.11 events at 173 GeV to 0.04 ±0.01 events at the kinematic limit. The 
resulting upper limits on cr(e + e _ — * XY) • BR(X — * Y7) for My ~ as a function of Mx range from 
0.46 pb to 0.075 pb, as shown in Figure |l0| 



4.2 Acoplanar-Photons 

The acoplanar-photons selection applied to the 130, 136 and 183 GeV data samples yields 2, 2 and 10 
events, respectively, in good agreement with the KORALZ predictions of 1.02 ± 0.02, 1.22 ± 0.02 and 
9.14 ±0.09 events for the Standard Model e + e~ — > 1/1/77(7) contribution. The expected contributions 
from other Standard Model processes and from cosmic ray and beam related backgrounds are small: 
less than 0.05, 0.05 and 0.1 events, respectively. The numbers of events expected and observed at 
the three centre-of-mass energies are summarized in Table ||. Also shown are the selection efficiencies 
for e + e~ — > 1/1/77(7) events, within the kinematic acceptance of the selection, and the corresponding 
cross-section measurements at each centre-of-mass energy. The OPAL measurements of the cross- 
sections at y/s = 130, 136, 161, 172 and 183 GeV are summarized in Table |4|. The measurements at 
130 and 136 GeV are weighted averages of the results obtained from the 1997 and 1995 data samples. 
The latter results, as well as the results at y/s = 161 and 172 GeV, have been taken from our previous 
publication [|J and are corrected for the different definition of the kinematic acceptance. 

Systematic errors arising from uncertainties on the electromagnetic calorimeter energy scale and 
resolution, the simulation of the detector material and consequent photon conversion probabilities, the 
integrated luminosity measurement and the detector occupancy estimate have been considered, and a 
relative systematic error of 8% is assigned to the cross-section measurements. This comes dominantly 
from uncertainty on the energy scale for low-energy photons and from comparison of different event 
generators |Q. 

The kinematic properties of the selected events in the combined 1997 data sample are displayed in 



Figures 11 and 12. They are compared with the predicted distributions for e + e~ — ► 1/1/77(7) obtained 
using the KORALZ generator normalized to the integrated luminosity of the data. In each case plot 
(a) shows the recoil mass distribution of the selected acoplanar-photon pairs. These are peaked near 
the mass of the Z° as expected for contributions from e + e~ — > 1/1/77(7). The resolution of the recoil 
mass is typically 2-4 GeV for M reco ii ~ Mz . Plot (b) shows the distribution of the scaled energy of the 
least-energetic photon. Plot (c) shows the 77 invariant mass distribution for which the mass resolution 
is typically 0.6-1.4 GeV. Plot (d) shows the distribution in scaled transverse momentum of the selected 
two-photon system. The measured kinematic properties of the events are given in Table ||. There were 
no selected events with three photons, compared to an expectation from KORALZ of 0.52±0.02 events. 

In data taken in 1995 at = 130 and 136 GeV, we observed 8 events compared to 1.6 ± 0.1 
expected from KORALZ |ffl|. However, the kinematic properties of these events, in particular the 
recoil mass distribution, agreed reasonably well with expectations. It was demonstrated that plausible 
detector or beam related backgrounds do not contribute to this excess. Analysis of the 1996 OPAL 
data samples taken at y/s = 161 and 172 GeV selected a number of events which was consistent 
with the expected Standard Model background contributions. In the data taken in 1997 at \fs = 
130 and 136 GeV, we select a total of 4 events where 2.24 ± 0.03 are expected. The 1997 data are 
therefore consistent with the KORALZ expectation. However, the OPAL data samples at y/s = 130 
and 136 GeV continue to favour an excess of events over the expectation from KORALZ. Results on 



14 



this topology have also been reported by the ALEPH collaboration based on the analysis of 5.8 pb -1 
of data taken at y/s = 130 and 136 GeV in 1995. That analysis has a kinematic acceptance similar 
to the one used by OPAL. No events were observed compared to a Monte Carlo expectation of two 
events [0211. 



As discussed in our previous publication, the status of event generators and analytical calculations 
of the Standard Model process e + e~ — > 1/1/77(7) is not yet satisfactory. We anticipate that on-going 
theoretical work by several authors will result in increased understanding of the actual precision of 
current approaches, and lead to improved approaches. In particular, a new and more complete calcula- 



tion has recently appeared [31]. For now, however, with contributions from higher order processes not 
demonstrably under control, we do not know what theoretical uncertainty to assign to the KORALZ 
prediction. 

In conclusion, the observed data excess over the KORALZ prediction at y/s = 130 and 136 GeV 
has not been resolved satisfactorily. Analysis of all data collected at these centre-of-mass energies 
by the other LEP experiments could help to resolve whether the observed effect is real, rather than 
simply a statistical fluctuation or a deficiency in the calculation of the Standard Model prediction. 

In the data taken at y/s = 183 GeV, the agreement with expectations from KORALZ is rather 
good. The only point of note is the selection of two events with a rather high invariant mass for the 
77 pair. As seen in Figure [T^ c these events populate the high-mass tail of the distribution expected 
from e + e~ — > 1/1/77(7). These events are both characterized by almost back-to-back photons in the 
polar angle region above | cos 6 \ > 0.9. 



4.2.1 Search for e+e - — > XX, X — »■ Y7 ; General case: My > 



Selected events are classified as consistent with a given value of Mx and My if the energy of each of 
the photons falls within the region kinematically accessible to photons from the process e + e~ — ► XX, 
X — * Y7, including resolution effects. The selection efficiencies at each generated grid point for the 
e + e~ — > XX, X — ► Y7 Monte Carlo events at y/s = 183 GeV are shown in Table |6|. These values 
include the efficiency of the kinematic consistency requirement which is higher than 95% at each 
generated point in the (Mx, My) plane. 

Figure 13 shows the 95% CL exclusion regions for cj(e+e- XX) • BR 2 (X Y7). The limits 
vary from 0.08 pb to 0.37 pb for M x > 45 GeV and M x - My > 5 GeV. In the region 2.5 GeV 
< Mx — My < 5.0 GeV, the efficiency falls off rapidly (see Table ^). However, even accounting 
for increased uncertainty on the efficiency, the limits in that region are better than about 1 pb for 
all (Mx,My). Because of the uncertainties in the modelling of the Standard Model background, as 
discussed earlier and in these limits and the limits presented below for this topology have been 
calculated without taking into account the background estimate. Events from e + e~ — * 1/1/77(7) are 
typically characterized by a high-energy photon from the radiative return to the Z° and a second 
lower energy photon. The kinematic consistency requirement is such that the two photons must 
have energies within the same (kinematically accessible) region. As Mx and My increase, the allowed 
range of energy for the photons narrows and fewer 1/1/77(7) events will be accepted. For the 10 selected 
events at y/s = 183 GeV, the distribution of the number of events consistent with a given mass point 
(Mx,My) is consistent with the expectation from e + e~ — * 2/1/77(7) Monte Carlo. 

Systematic errors are due primarily to limited Monte Carlo statistics at the generated (Mx,My) 
points and the uncertainty on the efficiency parametrization across the (Mx, My) plane. The combined 
relative uncertainty on the efficiency varies from about (3-6)% across the plane (for Mx — My > 5 
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GeV). All systematic uncertainties are accounted for in the manner advocated in reference |3(J]. This 
also applies to the limits for the My ~ case, presented in the next section. 



4.2.2 Search for e+e - — > XX, X — > Y7 ; Special case: My « 



For the special case of My ~ the kinematic consistency requirements differ from those used for 



the general case. One can calculate [ 10 1 the maximum mass, M^ ax , which is consistent with the 
measured three-momenta of the two photons, assuming a massless Y. A cut on M^ 1 ax provides further 
suppression of the 1/1/77(7) background while retaining high efficiency for the signal hypothesis. This 
is discussed in more detail in reference [||. We require that the maximum kinematically allowed mass 
be greater than Mx — 5 GeV, which retains 97 ± 2% relative efficiency for signal at all values of 
Mx while suppressing much of the remaining 1/1/77(7) background. Figure [lj|a shows the expected 
M™ ax distribution for signal Monte Carlo events with Mx = 80 GeV and for e + e~ — * 1/1/77(7) Monte 
Carlo events. Also shown is the distribution of the selected data events. In addition, we require that 
the recoil mass be less than 80 GeV, which approximately maximizes the expected sensitivity of the 
analysis for all Mx- This cut retains more than 50% of the signal efficiency, for all values of Mx, and 
dramatically reduces the residual 1/1/77(7) background. This is most true in the region of low Mx and 
remains valid up to the kinematic limit. A recoil mass cut is not applied in the massive Y case since 
this would lead to a large loss of efficiency in certain regions of the (Mx, My) plane. For the My ~ 
case, the efficiencies calculated from Monte Carlo events generated at 183 GeV are shown in Table |7| 
after application of the event selection criteria and the cut on M^ 12 ™, and then after the additional 
requirement of M reco ii < 80 GeV. Also shown in Figure [L4| are the M reco ii distributions for selected 
events with (b) no cut on M^ ax and (c) M^ ax > 75 GeV (for consistency with M x = 80 GeV). No 
event survives the recoil mass cut; the expected number of Standard Model events is 0.66 ± 0.03. The 
expected number consistent with Mx > 45 GeV is 0.36 ±0.02 decreasing to 0.07±0.01 expected events 
consistent with Mx > 90 GeV as shown in Table [7[ 

Based on the efficiencies and the number of selected events, we calculate a 95% CL upper limit 
on cr(e + e~ — ► XX) • BR 2 (X — > Y7) for My ~ as a function of Mx- This is shown as the solid line 
in Figure The limit ranges from 0.094 to 0.14 pb for all Mx from 45 GeV to the kinematic limit. 
Also shown as a dashed line is the expected limit, defined in section |4.1.1| . These limits can be used to 
set model-dependent limits on the mass of the lightest neutralino in supersymmetric models in which 
the NLSP is the lightest neutralino and the LSP is a light gravitino (X = x?,Y = G). Shown in 



Figure 15 as a dotted line is the (Born-level) cross-section prediction from a specific light gravitino 
LSP model ]l2| in which M{en) = 1.35M(x?), M(e£) = 2M(cr) and the neutralino composition is 
purely gaugino (bino). Within the framework of this model, Xi masses between 45 and 83 GeV are 
excluded at 95% CL. The expected number of 2/1/77(7) background events consistent with Mx > 83 
GeV is 0.12 ±0.01. 



As described in section 2, the efficiencies over the full angular range have been calculated using 
isotropic angular distributions for production and decay of X. The validity of this model has been 
examined based on the angular distributions calculated for photino pair production in ||] . For models 
proposed in ||, the production angular distributions are more central and so this procedure is con- 
servative. For a 1 + cos 2 6 production angular distribution expected for t-channel exchange of a very 
heavy particle according to ||, the relative efficiency reduction would be less than 2% at all points in 
the (Mx, My) plane. 
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5 Conclusions 



We have searched for photonic events with large missing energy in two topologies in data taken with 
the OPAL detector at LEP, at centre-of-mass energies of 130, 136 and 183 GeV. 

In the single-photon selection, which requires at least one photon with xt > 0.05 in the region 
15° < 6 < 165° (| cos 6 | < 0.966), a total of 21, 39 and 191 events are observed in the data for 

= 130, 136 and 183 GeV, respectively. These numbers are in agreement with the expectations of 
the KORALZ Monte Carlo generator for the Standard Model process e + e~ — ► 1/1/7(7). The expected 
background is small. We derive upper limits on the cross-section times branching ratio for the process 
e + e~ — > XY, X — * Y7 for the general case of massive X and Y. The limits vary from 0.075 to 0.80 
pb in the region of interest in the (Mx, My) plane and include the special case of My ~ 0, where the 
limit varies between 0.075 and 0.46 pb for the Mx mass range from M% to 183 GeV. 

The acoplanar-photons selection requires at least two photons with scaled energy x 7 > 0.05 within 
the polar angle region 15° < 6 < 165° or at least two photons with energy > 1.75 GeV with 
one satisfying |cos#| < 0.8 and the other satisfying 15° < 9 < 165°. In each case, the requirement 
pJ 7 /-Ebeam > 0.05 is also applied. A total of 2, 2 and 10 events are selected from the data samples at i/s 
= 130, 136 and 183 GeV, respectively. The KORALZ predictions for the contributions from e + e~ — ► 
1/1/77(7) are > respectively, 1.02, 1.22 and 9.14 events. The number of events observed in the 1997 
data samples and their kinematic distributions are consistent with expectations for e + e~ — * 1/1/77(7). 
We derive 95% CL upper limits on cr(e + e- -» XX) • BR 2 (X -► Y7) ranging from 0.08 to 0.37 pb for 
the general case of massive X and Y. For the special case of My ~ 0, the 95% CL upper limits on 
a(e + e~ — ► XX) • BR 2 (X — > Y7) range from 0.094 to 0.14 pb. Due to the uncertainties in the current 
modelling of the Standard Model process e + e~ — > 1/7/77(7), an limits from the acoplanar-photons 
analysis were calculated without taking into account the background estimate. 

For the single-photon and acoplanar-photons search topologies, the general case of massive X and 
Y is relevant to supersymmetric models in which X = X2 an d Y = Xi> with \\ ~^ X?7 an d Xi stable. 
The special case of My ~ is of particular interest for single and pair production of excited neutrinos 
and for supersymmetric models in which the LSP is a light gravitino and the NLSP is Xi which decays 
as Xi ~^ G7. For the latter scenario, the results of the acoplanar-photons search are used to place 
model-dependent lower limits on the Xi mass. A specific light gravitino LSP model [^] is excluded 
for the case of promptly decaying neutralinos with masses between 45 and 83 GeV. 
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Background process 


\/s=IoU GeV 


ys=13b GeV 


ys=183 GeV 


4- — -t- — / ^ 

e^e — > /ir // (7) 
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e T e — > ( T t ur/yy) 


no _i_ m 
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.Do ± .01 
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.lz ± .Uo 


e e / / / V / / 


01 4- 01 


01 4- 01 


no _l m 
.UO ZIZ .u± 


e + e~ — > 77(7) 


<.02 


<.03 


.03 ± .02 


Total physics bkgd. 


.07 ± .02 


.09 ± .02 


.38 ± .05 


Other bkgd. 


.02 ± .02 


.02 ± .02 


.39 ± .14 


Total background 


.09 ± .03 


.11 ± .03 


.77 ± .15 



Table 1: Numbers of events expected from various background processes contributing to the single- 
photon event sample for the three centre-of-mass energies. The Standard Model background contri- 
butions are given by process. Also shown are the expected other backgrounds coming from cosmic 
rays and beam related sources. The errors shown are statistical and the upper limits are at 68% CL. 



v^(GeV) 


Nobs 


N^ 7 ( 7 ) 


N bkg 


£i/I7 7 ( 7 )(%) 


CTmcas (Pb) 


130 


21 


25.8 ± 0.1 


0.09 ± 0.03 


77.0 ± 0.4 


11.6 ± 2.5 


136 


39 


31.2 ± 0.2 


0.11 ± 0.03 


77.5 ± 0.4 


14.9 ± 2.4 


183 


191 


201.3 ± 0.7 


0.77 ± 0.15 


74.2 ± 0.3 


4.71 ± 0.34 



Table 2: For each centre-of-mass energy, the table shows the number of events from the single- 
photon selection observed in the 1997 OPAL data, the number expected based on the KORALZ 
e + e~ — ► 1/177(7) event generator and the number of events expected from backgrounds. Also shown 
are the efficiencies for e + e~ — > 1/1/7(7) within the kinematic acceptance of the single-photon selection 
(defined in section 3) and the background-subtracted measured cross-sections within the kinematic 
acceptance. The errors shown are statistical. 



v^(GeV) 


Nobs 


N^j/ 77 ( 7 ) 


^i/^ 77 ( 7 ) (%) 


<rS (7, (pb) 


130 


2 


1.02 ± 0.02 


69.3 ± 1.0 


1.23 ±0.87 


136 


2 


1.22 ± 0.02 


69.1 ±0.7 


0.86 ±0.61 


183 


10 


9.14 ± 0.09 


67.9 ± 0.4 


0.27 ±0.09 



Table 3: For each centre-of-mass energy, the table shows the number of selected acoplanar-photons 
events in the 1997 OPAL data and the number expected based on the KORALZ e + e~ — ► 1/1/77(7) 
event generator. Also shown are the efficiencies for e + e~ — > 1/1/77(7) events within the kinematic 
acceptance of the acoplanar-photons selection (defined in section 3) and the corresponding cross- 
section measurements. 



v^(GeV) 


C7aZ (7) (pb) 


1/1/77(7) (rth \ 
KORALZ vP D J 


130 


1.49 ±0.68 


0.626 ±0.010 


136 


1.23 ±0.56 


0.526 ± 0.008 


161 


0.16 ±0.16 


0.330 ± 0.018 


172 


0.32 ±0.23 


0.303 ± 0.017 


183 


0.27 ±0.09 


0.247 ± 0.002 



Table 4: The measured cross-section for the process e + e~ — > 2/1/77(7), within the kinematic acceptance 
defined in section 3, for different centre-of-mass energies. For = 130 and 136 GeV the measurements 
are the weighted average of the results obtained from the 1997 data and the results obtained from the 
1995 data. The latter results, as well as those at y/s = 161 and 172 GeV, are taken from our previous 
publication [1] and have been corrected for the different definition of the kinematic acceptance. The 
final column shows the cross-section predictions from KORALZ. The quoted errors are statistical. 
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Xl 
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COS#i 


COS#2 


01 
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M 77 
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0.059 
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5.637 


5.006 
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29.8 


4.0 


0.218 


130.0 


0.263 


0.266 


0.296 


0.956 


3.574 


2.162 


91.4 


54.3 


20.0 


0.275 


135.9 


0.432 


0.063 


-0.801 


0.928 


3.037 


5.901 


99.1 


7.3 


22.2 


0.236 


136.0 


0.367 


0.234 


0.482 


0.657 


2.372 


0.847 


88.9 


57.2 


22.7 


0.373 


182.8 


0.484 


0.221 


-0.958 


0.125 


2.138 


1.318 


107.4 


74.6 


40.6 


0.329 


182.7 


0.674 


0.061 


0.867 


0.451 


4.367 


4.788 


94.7 


43.2 


11.8 


0.387 


182.7 


0.519 


0.518 


0.957 


-0.965 


4.848 


4.556 


84.1 


91.3 


91.2 


0.282 


182.7 


0.726 


0.096 


0.706 


-0.737 


5.076 


2.587 


90.3 


16.4 


47.1 


0.464 


182.7 


0.761 


0.026 


0.293 


-0.883 


4.921 


0.762 


87.2 


16.2 


22.4 


0.721 
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0.059 
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18.0 


20.3 
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182.6 


0.709 


0.067 


0.551 


0.844 


0.985 


1.005 


87.0 


45.5 


8.3 


0.627 


182.7 


0.541 


0.416 


-0.915 


0.901 


1.490 


6.115 


91.4 


82.7 


83.1 


0.271 


182.7 


0.541 


0.156 


-0.841 


0.903 


3.050 


5.363 


113.2 


25.3 


51.9 


0.252 



Table 5: Kinematic properties of selected acoplanar-photons events. Energies and masses are in GeV. 
Angles are in radians. The quantity M^ ax is defined in section 4.2.2. 



M x 


M Y =0 


M Y = M x /2 


M Y = M x - 10 


M Y = M x - 5 


M Y = M x - 2.5 


90 


71.8 ± 1.3 


71.0 ± 1.4 


65.5 ± 1.4 


39.2 ± 1.5 


6.8 ±0.8 


85 


72.5 ±1.3 


71.8 ± 1.3 


64.4 ± 1.4 


40.0 ± 1.5 


7.6 ±0.8 


80 


71.9 ± 1.3 


71.2 ± 1.4 


65.7 ± 1.4 


41.4 ± 1.5 


5.9 ±0.7 


70 


71.5 ± 1.4 


70.0 ± 1.4 


63.9 ± 1.4 


46.0 ± 1.5 


6.6 ±0.7 


60 


73.1 ± 1.3 


74.2 ± 1.3 


64.3 ±1.4 


43.5 ± 1.5 


7.7 ±0.8 


50 


73.7 ± 1.3 


72.2 ± 1.3 


64.3 ± 1.4 


43.7 ± 1.5 


9.4 ±0.9 



Table 6: Acoplanar-photons selection efficiencies (%) for the process e + e~ — ► XX, X — > Y7 at y/s = 183 
GeV for various Mx and My (in GeV) . These values include the efficiency of the kinematic consistency 
cuts. The efficiencies for the generated points at My = 20 and My = Mx — 15 are not shown. 





Selection efficiency with 


Selection efficiency with 


-^"1/^77(7) 


M x 


Mg™ > M x - 5 GeV 


M recoi i < 80 GeV 




90 


71.3 ± 1.4 


57.5 ± 1.5 


0.07 ±0.01 


85 


71.6 ± 1.4 


53.8 ±1.5 


0.10 ±0.01 


80 


70.0 ± 1.4 


47.1 ± 1.5 


0.13 ±0.01 


70 


69.2 ± 1.4 


40.7 ± 1.5 


0.18 ±0.01 


60 


70.4 ± 1.4 


42.7 ± 1.5 


0.25 ±0.02 


50 


70.4 ± 1.4 


37.2 ± 1.5 


0.31 ±0.02 



Table 7: Acoplanar-photons event selection efficiency (%), as a function of mass, for the process 
e + e~ — ► XX, X — > Y7, for My ~ at yfs = 183 GeV. The first column shows the efficiency of the 
selection described in section 3.2, after the cut on M^ ax - The second column shows the efficiency (%) 
after the additional requirement that M reco ii < 80 GeV. The last column shows the expected number 
of events from the process e + e~ — > 1/1777(7) (KORALZ). The errors are statistical. 
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Figure 1: For single-photon candidates in the data, a) shows the difference between the measured 
time-of-flight in the TOF and that expected for a photon from the interaction point for events passing 
all cuts or failing only the TOF timing cut; b) the same plot but magnifying the expected region for 
good events; c) shows the difference between the measured time-of-flight in the TE scintillators and 
that expected for a photon from the interaction point for events passing all cuts or failing only the 
TE timing cut. 
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Figure 2: a) The distribution for the cluster shape variable S, for the KORALZ e + e~ — > ^7(7) 
Monte Carlo (histogram) and for data (points). Events must have the primary photon candidate 
within I cos# | < 0.72 and passing all cuts or failing only the special background vetoes. The Monte 
Carlo sample is normalized to the integrated luminosity of the data, b) The transverse momentum 
sum of hits measured in the FD and SW calorimeters is plotted against the difference in azimuth 
between the forward momentum sum and the direction opposite the measured photonic system. The 
events plotted are those passing all cuts or failing only this cut. The solid circles are from the 
e + e~ —> e + e~7 Monte Carlo while the open triangles are from the e + e~ — > 1/177(7) Monte Carlo. Most 
e + e~ — ► 1/1/7(7) events have no significant forward energy and are not shown in the plot. The number 
of events correspond to ten times the integrated luminosity of the data. 
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Figure 3: The measured value of <r(e + e - — ► 7(7) + invisible particle(s)), within the kinematic 
acceptance of the single-photon selection, as a function of yfs. The data points with error bars are 
OPAL measurements at yfs = 130, 136, 161, 172 and 183 GeV. The curve is the prediction for the 
Standard Model process e + e~ — ► ^7(7) from the KORALZ generator. The data points at 130 and 
136 GeV represent the weighted means of cross-section measurements obtained from the 1995 and 
1997 data samples. The cross-section measurements from the previous data sets at y/s = 130, 136, 
161 and 172 GeV have been corrected for the difference in kinematic acceptance from the present 
analysis. 
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Figure 4: a) Distribution of x 7 vs cos 9 for the most energetic photon in the single-photon selec- 
tion at yfs = 183 GeV. The fine points are the KORALZ e + e~ — > 1^7(7) Monte Carlo (arbitrary 
normalization) and the solid triangles are the data, b) The cos# distribution for the most energetic 
photon in the single-photon selection at ^/s = 183 GeV. The points with error bars are the data and 
the histogram is the expectation from the KORALZ e + e~ — > uu^{^) Monte Carlo normalized to the 
integrated luminosity of the data. 
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Figure 5: The recoil mass distribution for events passing the single-photon selection for the y/s = 
130, 136 and 183 GeV data samples. The points with error bars are the data and the histograms 
are the expectations from the KORALZ e + e~ — > 1/177(7) Monte Carlo normalized to the integrated 
luminosity of the data. 
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Figure 6: Single-photon selection efficiency contours (in %) for ^/s = 183 GeV as a function of Mx 
and My. Lines are drawn around the boundaries defined by Mx + My = 183 GeV, Mx = My, and 
Mx + My = Mz, and to display the boundary between the small and large My regions. The shaded 
region indicates where the radiative return cut is applied. The region in which Mx — My is small is 
not considered, as explained in the text. 
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Figure 7: Number of single-photon candidate events in the y/s = 183 GeV data sample consistent 
with each set of mass values (Mx, My) for the process e + e _ — * XY, X — ► Y7. The boundaries and 
delineated regions are as defined for Figure |6[ 
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Figure 8: Expected number of single-photon events from the process e + e~ — > 1/1/7(7) at y / s=183 GeV 
which are consistent with each set of mass values (Mx, My) for the process e + e~ — > XY, X —* Y7. 
This figure gives the expected Standard Model contribution to Figure |7[ KORALZ was used to model 
the e + e~ — ► 1/1/7(7) process. The boundaries and delineated regions are as defined for Figure ||. 



29 



OPAL 




M x (GeV) 

Figure 9: The 95% CL upper limit on fj(e+e" -> XY) • BR(X -» Y7) at yfa = 183 GeV as a function 
of Mx and My. The boundaries and delineated regions are as defined for Figure ||. 
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Figure 10: The 95% CL upper limits on <r(e + e- -► XY) -BR(X -► Y7) at ^fs = 183 GeV as a function 
of Mx, assuming My ~ 0. 
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Figure 11: Plots of kinematic quantities for the selected acoplanar-photons events from the combined 
\fs = 130 — 136 GeV data sample, a) Recoil mass distribution, b) Distribution of the scaled energy of 
the second photon (X2). c) Distribution of the invariant mass of the 77 system, d) Scaled transverse 
energy distribution for the 77 system. The data points with error bars represent the selected OPAL 
data events. In each case the histogram shows the expected contribution from e + e _ — ► ^1777(7) events, 
from KORALZ, normalized to the integrated luminosity of the data. 
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Figure 12: Plots of kinematic quantities for the selected acoplanar-photons events for y/s = 183 
GeV. a) Recoil mass distribution, b) Distribution of the scaled energy of the second photon (X2). c) 
Distribution of the invariant mass of the 77 system, d) Scaled transverse energy distribution for the 
77 system. The data points with error bars represent the selected OPAL data events. In each case the 
histogram shows the expected contribution from e + e~ — > 1/1/77(7) events, from KORALZ, normalized 
to the integrated luminosity of the data. 
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Figure 13: The shaded areas show 95% CL exclusion regions for cr(e+e- -»• XX) • BR 2 (X -► Y 7 ) at 
\fs = 183 GeV. No limit is set for mass-difference values Mx — My < 5 GeV, defined by the lower 
line above the shaded regions. The upper line is for Mx = My. 
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Figure 14: Properties of selected acoplanar-photons events at ^fs = 183 GeV. a) The M™ ax distri- 
bution, b) The recoil mass distribution prior to any Aff requirement, c) The same recoil mass 
distribution with a cut at M^ ax > 75 GeV (this cut is for consistency with an Mx of 80 GeV). In 
each plot the OPAL data are shown as points with error bars, the shaded histogram shows the ex- 
pected distribution for e + e~ — > uT7jj(j), normalized to the integrated luminosity of the data, and the 
unshaded histogram shows the expected distribution for the signal process e + e~ — ► XX, X — > Y7 for 
Mx = 80 GeV. The three distributions for signal Monte Carlo are normalized to the same (arbitrary) 
production cross-section. 
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Figure 15: 95% CL upper limit on er(e + e~ -» XX) • BR 2 (X -» Y7) for My ~ (solid line). Also 
shown is the expected limit (dashed line). The dotted line shows the cross-section prediction of a 
specific light gravitino LSP model (T^j. Within that model, Xi masses between 45 and 83 GeV are 
excluded at the 95% CL. These limits assume that particle X decays promptly. 



36 



